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I. INTRODUCTION 



One of the basic ingredients of the Standard Model is the Conserved Vector Current hypothesis (CVC), which 
postulates that the charged (isovector) weak vector current and the neutral isovector component of the electromagnetic 
(EM) current are members of the same isovector multiplet. Since the charged current spectral function is now 
measured rather accurately in r decay [1-3], CVC can be tested experimentally, provided, that is, one can extract the 
isovector spectral function from data on e + e~ — * hadrons [4-8]. In the absence of isospin breaking, this extraction is 
straightforward since, for example, for n-pion final states, a state with an even (odd) number of pions has even (odd) 
G-parity and hence can be produced only through the isovector (isoscalar) component of the hadronic EM current. 

Isospin breaking, however, complicates the extraction of both the isovector and isoscalar spectral functions. Before 
proceeding, it is useful to clarify our notation. We define the standard SU(3) f octet of vector currents by = q^/^^-q, 
where A a are the usual Gell-Mann matrices. The electromagnetic current is then J E M = J 3 + J^/y/H, while the scalar 
correlators, H ab (q 2 ) (where we will restrict our attention to a, b — 3, 8) are defined by 

i J d 4 x exp(igz) (0|T (J» J») |0) = - <?V) n Q V). (1) 

Defining the spectral functions, p ab (q 2 ), corresponding to the II ab (q 2 ) in the standard manner, p ab (q 2 ) = -^ImH ab (q 2 ), 
the EM spectral function is then 

P EM = P 33 + ^P 38 + Ip 88 - (2) 

p EM contains the isovector spectral function, p 33 , which is the isospin rotated version of the corresponding charged 
current isovector spectral function measured in r decays. Owing to the relation, 

a(e+e- — hadrons) = ^- p EM ( s ), (3) 

s 

however, the cross-section for e + e~ — > hadrons directly measures p EM , and not p 33 or p 88 . If isospin were not 
explicitly broken, this would present no problem since p 38 would necessarily vanish and, as noted above, one could in 
addition identify the states contributing to the isovector (33) and isoscalar (88) spectral functions by their G-parity. 
Near threshold it is known, from a study of the mixed isospin vector correlator (ab = 38) to two loops in chiral 
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perturbation theory (ChPT) [9], that isospin breaking effects in p EM are negligible. However, in the region of the p 
and oj resonances, isospin breaking is significant, as signalled by the interference dip in the cross-section for the tttt final 
state [10]. The conventional method [5] for making corrections for this observed isospin breaking, in order to extract 
the vector component of the EM spectral function, is to first parametrize the amplitude in terms of a sum of p and 
u> Breit-Wigner resonance forms (in general one includes also contributions associated with the higher p resonances) , 
and having fitted the parameters to the observed cross-sections, remove the uj contribution to the amplitude by hand. 
The squared modulus of the resulting modified amplitude is then used in place of the squared modulus of the original 
amplitude to identify that portion of the cross-section to be associated with the purely isovector (33) portion of the 
EM spectral function. 

The conventional proceedure just described for correcting for isospin breaking, however, does not, in fact, produce 
the desired 33 component of the vector spectral function. To understand why this is the case, let us first define, for 
the neutral vector mesons, V = p, u, <j>, the decay constants Fy via 

(0\j;\V(k)}=rn v F v e,(k) (4) 
where e t _ l (k) is the vector meson polarization vector, and a = 3, 8. Fp , Ft and F^ are non-zero in the isospin limit, 

(S) (3) (3) 

while Fp , Ft, and F^ vanish in the isospin limit and hence are proportional to the isospin breaking parameter 
5m = rrid — m u . In the presence of isospin breaking, all of the neutral vector mesons, in principle, mix with one 
another, so the physical states are of mixed isospin. If we consider the tttt final state mediated by the u> exchange 
then, to leading order in isospin breaking, this transition is indeed mediated by the isoscalar component of the EM 
current (the intermediate u> contribution generated by the isovector current is second order in isospin breaking, one 

factor from the coupling Ft , and one from the isospin violating u> — > tt + tt~ decay vertex), and hence should be 
removed if one wishes to extract only the isovector contributions. The remaining p exchange contributions are, for a 
similar reason, however, not purely isovector. The reason is that if one considers the p exchange contribution to the 
EM spectral function, the 38 component is first order in isospin breaking (being proportional to FpFp ). Thus, the 
p contribution to e + e~ — ► -k + tt~ necessarily contains a piece first order in isospin breaking, and associated with the 
38 part of the EM spectral function, which must be subtracted in order to isolate the purely isovector 33 component. 
A similar argument shows that there will also be a 38 contribution to the measured cross-section for e + e~ — > u> —* 3ir, 
which one would have to correct for in order to isolate the purely isoscalar 88 component of p EM . 

In what follows we will discuss how to perform the corrections associated with first order isospin breaking contri- 
butions to the vector meson decay constants. The paper is organized as follows. In section II we first discuss some 
general issues, which provide useful qualitative guidelines for the subsequent discussion. In addition we discuss an 
analogous case involving the neutral isoscalar and isovector axial vector currents, which example serves to illustrate 
the basic features we will meet in the vector current case of interest, but in a context where, unlike the vector case, 
the isospin breaking decay constants are already known, being fixed by a next-to-leading order analysis using Chiral 
Perturbation Theory (ChPT). In section III, we then show how, and with what accuracy, existing QCD sum rule 
analyses of the mixed isospin (38) vector current correlator can be used to extract the isospin breaking decay constants 

(8) {3} {3} 

F p Ff and F^> . In section IV, we use the results of this analysis to evaluate the corrections required to extract 
the pure isovector contribution associated with the p, and pure isoscalar contributions associated with the u> and 4> 
from the EM cross-section data, and comment on the effect such corrections would have, for example, on precision 
tests of CVC and the extraction of chiral low-energy constants (LEC's) via the inverse chiral sum rules method [11]. 



II. SOME GENERALITIES, AND AN ILLUSTRATIVE EXAMPLE 

In the Introduction we have explained the reason for the existence of previously neglected isospin breaking cor- 
rections in the extraction of the isoscalar and isovector spectral functions from e+e~ — ► hadrons data. We have, 
however, not yet demonstrated that such corrections can be expected to be numerically significant. As will be seen 
below, making numerical estimates for the size of these corrections is non-trivial, and we will be forced to rely on a 
QCD sum rule analysis of the mixed-isospin vector current correlator to make these estimates. Since it can be difficult 
to estimate, in a quantitative manner, the errors present in such an analysis, it is useful to consider any qualitative 
constraints, based on general principles, which one might use to judge the plausibility of the resulting solutions. We 
will discuss below the existence of such constraints based on the framework of effective chiral Lagrangians. We will 
also consider, as an illustrative example, an analogous case in which the vector currents of the problem at hand are 
replaced by axial vector currents. The advantage of this example is that the isospin breaking pseudoscalar decay 
constants (anologous to the isospin breaking vector meson decay constants to be determined by the QCD sum rule 
analysis below) can, in this case, be computed with good accuracy using the methods of ChPT. This allows us to show 
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explicitly, in a context where the numerical accuracy of the evaluation of the isospin breaking corrections is not open to 
question, that such isospin breaking corrections can play a significant role in the correct extraction of flavor-diagonal 
spectral contributions from the spectral functions of correlators involving mixed- flavor currents. Certain features of 
the relation of the relative signs and magnitudes of the corrections in the isovector and isoscalar cases, which recur 
in the vector channel, will also be exposed, again in a context where the accuracy of the numerical estimates is not 
open to question. From this example we will be able to unambiguously conclude that isospin breaking corrections of 
the type also present in the vector channel must be expected to be numerically significant, especially for observables 
related to differences of weighted integrals of the isovector and isoscalar spectral functions, for which there will be 
cancellation between flavor-conserving contributions, but coherence between the isospin breaking corrections. 

Let us turn then to the qualitative guidance offered, in the vector channel, by the framework of effective chiral 
Lagrangians. Certain qualitative features of mixing in the vector meson sector follow immediately from the properties 
of an effective chiral Lagrangian such as would be obtained by the Callan, Coleman, Wess, Zumino [12,13] construction 
[14,15]. Note that the lowest order term in quark masses and derivatives which produces isospin mixing in the vector 
meson propagator matrix involves a single power of the quark mass matrix (where, as per the usual chiral counting, 
external momenta are counted as 0(q) and quark masses as 0{q 2 )). Moreover, at leading order, the isospin violating 
decay constants associated with the isospin pure states vanish. As a result, at this order, the transformation from the 
original isospin pure basis to the physical, mixed isospin basis, is a rotation, and the isospin breaking decay constants 
of the physical particles result purely from the mixing in the physical states. If we concentrate on p-u> mixing, 

(8) (3) 

this would mean that Ff> and F^j were equal in magnitude and opposite in sign. If we consider effective isospin 
breaking operators higher order in the quark mass/derivative expansion, new effects come into play. First, one finds 
non- vanishing "direct" isospin violating couplings of the external vector currents to the isospin pure states from terms 
involving both derivatives and one power of the quark mass matrix. Second, terms involving both derivatives and one 
power of the quark mass matrix can produce off-diagonal mixing elements in the wavefunction renormalization matrix, 
a consequence of which is that the transformation from the isospin pure to physical basis is no longer a rotation, but 
rather the product of a symmetric matrix and a rotation. Third, terms with two powers of the mass matrix will 
produce modifications in the momentum-independent mixing terms in the vector meson propagator matrix. One 
would expect the effect of such higher order terms to be manifest in deviations from leading order relations such as 
Fp 8 ^ = — F^ 3 \ The size of such deviations should be typical of next-to-leading order corrections in SU(3)l x SU(3)r, 
and hence might be as large as <~ 30%. 

We next turn to the illustrative example mentioned above, in which we replace the vector mesons and vector 
currents with pseudoscalar mesons and axial vector currents. Not only can several of the qualitative points just made 
be clearly illustrated in this case, but the actual numerical values of the relevant isospin breaking corrections can, for 
this example, be calculated to good accuracy using the techniques of ChPT, since all of the relevant decay constants 
are known at next-to-leading order in the chiral expansion. This is in contrast to the case of the vector current spectral 
functions, where we will be forced to rely on a QCD sum rule analysis to obtain the isospin breaking decay constants. 
Let us, therefore, define the axial current combination 

= ^ + ~^A S ^ (5) 

where A 3 / 8 are the 3, 8 members of the usual axial current octet. We then define the pseudoscalar decay constants, 
fp, for P = 7T°, i] and a = 3, 8, via 

(Q\Al\P{k))=if a P k». (6) 

At leading (second) order in the chiral counting the physical 7r°, rj basis is a pure rotation of the isospin pure octet 
basis 7r 3 , 7r 8 

7T° = 7T 3 + <9 , 7T 8 T, = 7T 8 - <9 7T 3 (7) 

where 6 a = \/3(md — m u )/A{m s — fa), with m = (m^ + m„)/2, is the leading order mixing angle, and the isospin 
breaking decay constants are produced purely due to the "wrong-isospin" admixture present in the physical tt°, rj 
states, 

fW=0 o F, fW = -6 F (8) 

(3^ (8) 

where F is a second order LEC, identical to both = f n , the physical pion decay constant, and = /,,, the 
physical r] decay constant, at this order in the chiral expansion. When one considers the full effective Lagrangian at 
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next-to-leading order, however, one finds that (1) there is an infinite renormalization required to regularize the mixing 
of 7r 3 , 7r 8 even if one works with the versions of the fields renormalized in the isospin limit; (2) there is indeed an 
off-diagonal element produced in the wavcfunction renormalization matrix and (3) there are indeed "direct" isospin 
breaking meson-current vertices. (For a detailed exposition of the above features beyond leading order see Ref. [16].) 
The effect of all these features is to produce the following results for the isospin breaking decay constants [17], written 
in terms of the next-to-leading order expressions for the tt° and r\ decay constants, which differ at this order in the 
chiral expansion, 



4 8) = e,U 

4 3) = 



(9) 



In Eq. (9), the quantities €12 differ by terms which are next-to-leading order in the chiral expansion (the complete 
expressions can be found in Ref. [17]) and, using the observed experimental ratio of /k/ fn, the next-to- leading order 
expressions for the isospin conserving decay constants imply f v / f n ~ 1.3 [17]. 

We are now in a position to consider the analogue of the vector current case of interest. To this end we imagine 
that we would like to obtain the ir° and r\ contributions to the isovector and isoscalar axial spectral functions of the 
scalar correlators, n 33 (<7 2 ) and I±f 8 (g 2 ), defined by 

(10) 



TL% = 1 J d 4 xe^(0\T(A;(x)At(0))\0) = Ilf (q 2 )q^q v + Hf(q 2 ) (q 2 g, v - q^q v ) . 

(We concentrate on the scalar correlators, H\ b , since it is these correlators which contain the pole contributions 
analogous to those of the vector mesons in the vector current correlators.) It is straightforward, in this case, to simply 
evaluate the contributions to the spectral functions, pf 3 . Keeping only terms up to first order in isospin breaking and 
to next-to-leading order in the chiral expansion, one finds that 



p?V) = 


r f(3)" 

J IT 


2 S(q 2 - ml) 


pf{q 2 ) = fl 3) fi 8) Hi 2 rnl) + /f/f S (q> m 


pfV) = 


/f" 


S(q 2 -m 2 ). 



(11) 

Let us imagine, however, that the way we had to go about extracting these contributions was by analyzing the 
"experimental" spectral function Pi(q 2 ) of the scalar correlator, U^, appearing in the analogue, IT^, of the correlator 
of the product of two EM currents, i.e., 



K> = i J d 4 xe^(0\T(A^(x)A v (0))\0) ^nf( q 2 )q^ + Ii^{q 2 ) (q 2 g^ - q^) . 



(12) 



The analogue of the standard vector current extraction of the isovector component would then consist of identifying 
the 7T° pole contribution with the isovector p\ 3 component of pf and the 77 pole term with the isoscalar p\ 8 component 
thereof. We can see, however, that this identification is incorrect. Indeed, the n° contribution to pf, to first order in 
isospin breaking, is straightforwardly found to be, not p 33 , but 



[P?l = P? + ^= 3 [Pfl - + ^ 3) 4 8) J % 2 - ml). 

Similarly, the -q pole contributions to pf, to first order in isospin breaking, is not ^pf 8 but 

[pi\ = \pf + ^[pl\={^ 



(13) 



f(8) 
J n 



2 + ^)f^\s(q 2 -ml). 



(14) 



Thus, to "extract" the isovector and isoscalar spectral functions from the "experimental" spectral function one would 
actually have to know the isospin violating decay constants and then form the combinations 



Pf =[ P tl-^fi 3) f^^-ml) 
pf = 3[pt] v -2V3tf \f^ S(q 2 -m%). 



(15) 
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Note that, because of the structure of the original current, the size of the isospin breaking correction in the isoscalar 
case is naturally larger by a factor of 3 than that in the isovector case. If we work with the values of ei t 2 from Ref. [17] 
(which correspond to the value of the isospin breaking mass ratio r = (mj — m u )/(md + m u ) used in the sum rule 
analysis we will employ below) 

ei = 1.37 x 1(T 2 , e 2 = l.llxl(r 2 (16) 

we find that that the "experimental" pion pole contribution to p^ is a factor of 1 + 2ei/\/3 = 1.016 larger than 
the true isovector spectral function, and that the nominal "experimental" isoscalar spectral function, obtained by 
taking the rj pole contribution and multiplying by 3, is smaller than the true isoscalar spectral function by a factor of 
1 — 2\/3e2 = 0.961. Thus, to extract the true isovector and isoscalar spectral functions, one would have to multiply the 
nominal "experimental" ones by 0.984 and 1.040, respectively. Note that the corrections go in the opposite direction 
for the two cases and that the magnitude of the correction is significantly larger in the isoscalar case. The reason for 
the latter feature of the results has already been explained. The reason for the former is that the leading order result 
that the ir° and r\ contributions to the 38 spectral function differ in sign, but not in magnitude, is still approximately 
satisfied by the next-to-leading order expressions, as expected based on the general arguments above. Note also that 
the effect of these: corrections can be greatly enhanced if one considers combinations which would vanish at leading 
order in the chiral counting. As an example, if we consider the integral over the 88-33 spectral function, the nominal 
value (without the above corrections) is 0.69/ 2 , while the actual value (after corrections) is 0.77/ 2 . The isospin 
breaking corrections thus represent a 12% increase for this quantity, much larger than one would guess based on the 
typical few percent size of familiar isospin breaking corrections. 

We will see in what follows that many of the features of the axial current example are recapitulated in the vector 
current case. In fact, because in the limit in which one considers the vector meson multiplet to be ideally mixed, but 
takes the decay constants to be otherwise determined by SU(3)p symmetry, the u EM decay constant is a factor of 3 
smaller than the p EM decay constant (rather than just the factor of \/3 difference between the "A" decay constants 
of the 7T° and the rj in the above example), the discrepancy between the size of the corrections in the isovector and 
isoscalar channels would be expected to be even greater in the vector case. We will see below that this expectation is 
indeed borne out. 



III. THE QCD SUM RULE EXTRACTION OF THE ISOSPIN BREAKING VECTOR MESON DECAY 

CONSTANTS 

Since the 3 and 8 components of the vector current octet occur in the standard model only in the combination 
J? M , it is not possible to directly determine the isospin breaking vector meson deay constants experimentally. One 
can, however, obtain indirect access to these quantities via a QCD sum rule analysis of the vector current correlator 

n 38 . 

The basic idea of such a QCD sum rule analysis [18-20] is straightfoward. From the behavior of U 38 (q 2 ) as q 2 — > oo 
in QCD, it is known that II 38 satisfies an unsubtracted dispersion relation. This dispersion relation allows one to relate 
the dispersion integral over hadronic spectral function p 38 to the value of II 38 at large spacelike value of q 2 = —Q 2 . 
The latter can be expressed in terms of vacuum condensates using the operator product expansion (OPE), while the 
hadronic spectral function depends on the isospin breaking and isospin conserving decay constants of the various 
vector mesons. The utility of this relation is greatly enhanced, as first noted in Ref. [18], if the original dispersion 
relation is Borel transformed since, in that case, the higher s portions of the transformed hadronic spectral integral are 
exponentially suppressed, while the contributions of higher dimensional operators on the OPE side are simultaneously 
factorially suppressed. In favorable circumstances one is then able to write the dispersion relation in a form in which 
the parameters of a small number of resonances dominate the hadronic side, and the contributions of a small number of 
vacuum condensates of low dimension operators (which condensates can be determined from other sum rule analyses) 
dominate the OPE side. One can then use the known values of the vacuum condensates to extract the (in our case 
unknown) resonance parameters. Note that it is far preferable, in the case of interest to us, to investigate the isospin 
breaking decay constants by analyzing an isospin breaking correlator, rather than by looking at the isospin breaking 
corrections to a correlator which also has an isospin conserving piece. This is because of the fact that, in the latter 
case, one would have to consider isospin breaking not only in the vector meson decay constants, but also in the 
continuum thresholds, which would be required in the spectral ansatze in order to model the higher s portions of the 
hadronic spectral function. Since that portion of the spectral model is a rather crude representation of the actual 
continuum, this would introduce potentially large, and difficult to control, uncertainties into such an analysis. 

Let us turn then to the sum rule analysis of the correlator II 38 (q 2 ). We will, in fact, use the results of existing 
analyses [21-24] of the related correlator, H pLU (q 2 ), defined by 
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n$(g) =ij d A x exp(uj • x)(0\TJP(x)J?(0)\0) 

= {q,q v -q 2 g^)W{q 2 ), (17) 

where = (u-f p u — d'y l _ l d)/2 and J" — (u'j^u + d'y^d)/^. This is possible because, truncating the OPE at 
0(m q ), 0(a s ), 0{oiem) and operators of dimension 6 (with either the vacuum saturation hypothesis, or the rescaled 
version thereof, for the four-quark condensate contributions), the ss portion of J® does not contribute to the OPE 
representation of II 38 (q 2 ). One then has 

[n 3 V)] OP£ = V3 [n-V)] OP£ , (is) 

where, after Borel transformation (indicated by the operator B) of the truncated OPE expression, one finds [21] 

b rwioPE = ^ h m 2 + Cl + ^ + ^ i , (19) 



with M the Borel mass, 



c 



OEM 

16tt 3 



ci = 0(m 2 ) ~ 

/ m„-m d (l+7)\ 

C2 = H 2T^ J (M>0 

224tt7 ._ l2 28tt ._ 
c 3 = gl a s K (qq} - -^-^ e K (ggr) , (20) 

where 7 = (dd) o/(uu)o — 1, (99)0 = ((S«)o + (dd)o)/2, a e and a s are the electromagnetic and strong coupling 
constants, respectively, and k is the parameter describing the deviation of the four-quark condensates from their 
vacuum saturation values. 

The hadronic spectral function, p 38 , is parametrized as 

P 38 ( S ) = ^ [f p 5 p (s) - Ms) + Uhis) + ■■■} (21) 

where 

Ms) = ^(- m ^ +m ^ (22) 

(which reduces to S(s — my) in the narrow width approximation), and we have followed the notation and sign 
conventions employed in the literature for the parameters fy describing the strengths of the various vector meson 
contributions to the spectral function in question. Note that, in the notation introduced above, 



fv - + F ( 3 M 8 ) 
4V3 



± F W F W (23) 



where the upper sign holds for V = p, (f>, • • •, and the lower sign for V = u>. 

A few comments are in order concerning the form of Eq. (21). The first concerns the presence of the 4> contribution, 
which was not included in the earliest sum rule analyses [21,22]. Recall that we expect, as a leading order result, that 
f p ~ f u (the absence of a minus sign in this relation is a consequence of the sign convention for the definition of f u ). 
This means that, since from far in the spacelikc region the p and u> masses appear essentially the same, one must 
expect significant cancellation between the p and w contributions to the original dispersion relation for large values of 
Q 2 . This is especially true in the narrow width approximation. Based on this observation, it was realized, in Ref. [23], 
that, although could be expected to be significantly smaller than f p or / w (of order 6-7% if one makes an estimate 
based on the Particle Data Group evaluation of the deviation of vector meson mixing from ideal mixing [23]), the 
contribution of the <j) to the actual sum rule need not be negligible. The sum rule was then re-analyzed, including 
the <f> contribution, though still in the narrow width approximation. The results supported the qualitative arguments 
regarding the importance of including the 4> contribution, and simultaneously cured a physically unpleasant feature of 
the earlier analyses, in which contributions from the p'-uj' region of the spectral function were as important, or more 
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important, than those from the p-u region in determining, through the original dispersion relation, the value of the 
correlator near q 2 = 0. It was, however, then subsequently pointed out [24] (IJL), again because of the high degree of 
cancellation between p and u> contributions to the dispersion integral in the far spacelike region, that the use of the 
narrow width approximation for the p might also be a rather poor one. This was borne out by the numerical analysis 
of Rcf. [24]. Of particular note is the fact that, introducing the p width into the spectral ansatz for the p contribution, 
one finds that f p and f w are lowered in magnitude by a factor of ~ 6, and by a factor of ~ 2 compared to the 
values extracted using the narrow width analysis. With the above understanding in mind, we will, therefore, in what 
follows, employ the spectral ansatz of Eq. (21) in the form arrived at in Ref. [24], i.e., with the ui and 4> treated in the 
narrow width approximation, but the p treated using the Breit-Wigner form given in Eq. (22). 

It should be stressed that, in using the results of Ref. [24], the present analysis relies only on the extraction of the 
hadronic spectral function for the mixed-isospin vector current correlator, accomplished in that reference using QCD 
sum rules, and not on the attempt by the authors of Ref. [24] to interpret this spectral function in terms of off-shell 
vector meson propagators. The latter interpretation (and that of Ref. [22]) is necessarily incorrect, since off-shell 
Green functions are well-known to be altered by redefinitions of the hadronic fields and, as such, are not capable of 
being related to physical objects such as the mixed-isospin vector current correlator (an earlier claim by the present 
author, contained in Ref. [23], that the rescaled versions of the vector fields represented a possible field choice for the 
vector mesons is also incorrect). The only questions relevant to the use of the spectral function solution of Ref. [24] in 
the present work are then (1) is the sum rule reliable for the scales at which it is employed and (2) is the original ansatz 
for the spectral function plausible in form. While it is not possible to provide a rigorous proof of the suitability of the 
analysis of Ref. [24], there is considerable indirect evidence in its favor. First, although a resonance-saturation ansatz 
involving s-dependent widths, but constant real parts of the resonance self-energies, does not rigorously implement 
the constraints of analyticity and unitarity, it is well-known from phenomenological studies of e+e~ — > hadrons that 
such ansatze, nonetheless, can provide fits of very high numerical accuracy to the experimentally measured spectral 
function (see Ref. [25] for a recent detailed discussion). This is also true of the resonance-saturation fits to the timclike 
pion form factor measured in hadronic r decays mediated by the isovector current (where, for example, the naive form 
of the resonance contributions involving s-dependence only in the widths, produces a fit of even slightly better quality 
than does an ansatz employing the Gounaris Sakurai form, which correctly implements analyticity and unitarity - see 
the results of Table 3 of Ref. [1]). Moreover, a recent analysis of the isovector (isospin conserving) vector correlator 
using a continuous family of finite energy sum rules as a means of implementing the dispersion relation implicit in 
QCD sum rules shows that fitting a naive form of the resonance saturation ansatz for the hadronic spectral function 
to the OPE representation (where the widths and decay constants of the resonances are free parameters determined 
by the fit to the OPE) reproduces the experimentally measured spectral function [1] to an accuracy of a few percent 
and, moreover, requires a phenomenological value of the p width within a few MeV of that obtained by direct fitting 
of the experimental data using the favored HLS model in Ref. [25] [26]. This indicates that, at least in the isospin 
conserving case, matching of an ansatz of the type employed in Ref. [24] to the OPE representation of the vector 
current correlator provides a good fit to the actual vector current spectral function. Note that the matching between 
the hadronic and OPE sides of the conventional QCD sum rule obtained in Rcf. [24] is also considerably better than 
that obtained in earlier analyses using the narrow width approximation for the p, further indicating the necessity of 
employing an ansatz of the IJL form. The only improvement in the IJL ansatz one might have hoped for was the use 
of a variable p width, to be fixed by the sum rule analysis. Since, however, in the isospin conserving case, the result 
of such an exercise is to essentially reproduce the width employed by IJL, it seems unlikely that the analysis would 
have been significantly altered by such an extension. 

In light of the discussion above, we should thus be able to extract the desired isospin breaking vector meson decay 
constants from the results of Ref. [24]. The results we use below, however, differ somewhat from those quoted in that 
reference, and the reasons for this difference must first be explained. The first point in need of clarification is related 
to the fact that there are two sets of results associated with the analysis using the physical p width in Ref. [24] . Of 
these two sets, only the one contained in the column labelled "physical widths" and "no constraint" in Table 1 of 
that reference should be employed. The reason for rejecting the other set (labelled "constrained") is that these results 
were obtained assuming II 38 (mj) could be extracted from the ui — > tttt contribution to the physical cross section. This 
assumption, however, as we have seen above, neglects the presence of additional II 38 contributions in the p exchange 
contribution to the cross-section, and hence is incorrect. The second modification we make in employing the results of 
Ref. [24] concerns the magnitude of the errors quoted on the extracted values of the parameters /y. It turns out that, 
in Ref. [24] an overly conservative error was assumed on the crucial input isospin breaking light quark mass ratio, 
r m = m u /md- The authors of Ref. [24] employed r m = 0.50 ± 0.25. The ratio, r m , however, is actually much better 
constrained than this from ChPT analyses [27]. As a consequence, the quoted errors in Ref. [24] are unnecessarily 
inflated. The authors of Ref. [24] have kindly provided unpublished results corresponding to the more realistic input 
fm = 0.54 ± 0.04 employed in earlier sum rule analyses of the correlator at hand. We will, therefore, employ, for 
determining the errors on the extraction of the desired isospin breaking vector meson decay constants, the results of 
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the "unconstrained" fit obtained using the modified input above, r m = 0.54 ± 0.04. The results, which correspond to 
a stable Borel regime 1.15 GeV < M < 2.45 GeV, are [28] 



f p = AVSF^F^ = 0.0030 ± 0.0012 GeV 2 
f u = -AVSF^F^ = 0.0025 ± 0.0009 GeV 2 
U = aV3F ( , 3) F ( , 8) = -0.0002 ± 0.0002 GeV 2 . 



(24) 



The contributions from the p' and w' are found to be very small, and cannot be reliably extracted from the sum rule 
in its current form. Note that the reliability of the results is supported by that fact that they display two features 
which correspond to our general expectations: first, that f p ~ f u and, second, that is of order 6-7% of f p and f w . 



We are now in a position to evaluate the isospin breaking decay constants Fp 
the results of Eq. (24), together with the relations 



(8) 



and Jtf>. 



To do so we employ 



t?EM _ P 3 , 



V3 



V ' 



(25) 



where the physical EM vector meson decay constants, Fy M , determined experimentally from the partial widths for 
the decays V — > e+e~, are 



Fj? M = 154 ±3.6 MeV 
F^ M = 45.9 ±0.8 MeV 
FF M = -79.1 ± 2.3 MeV. 



(26) 



The sign of the (f> decay constant in Eq. (26) has been chosen to be consistent with expectations from SU(3)f 
symmetry and ideal mixing. Note that the relative magnitudes in Eqs. (26) are roughly in line with the expectations 
of that limit (in which one would expect the p, lo and <p decay constants to be in the ratios 3:1: —y/2). It is then 
straightforward to solve for Fy and F v separately. One finds, for the isospin breaking decay constants, 

Ff ] = 2.8 ± 1.1 MeV 



F^> = -4.2 ± 1.5 MeV 
F] 3) = 0.21 ± 0.21 MeV, 



(27) 



where the quoted errors are totally dominated by the errors of the sum rule fit values of the products of decay 
constants. The values of the isospin conserving decay constants then follow immediately from Eqs. (25), (26) and 
(27). We will quote them in the form of ratios to the relevant experimental values, which form allows one to directly 
compute the additional corrections required to obtain the pure isovector and isoscalar spectral functions from those 
obtained conventionally, i.e., via the standard analysis described above. We find 



p(3) 

pEM 
P 



F 



(8) 



VSF L 



EM 



F 



(8) 



V3> 



EM 



- 1 = -0.011 ±0.0043 



1 = 0.091 ± 0.029 



1 = 0.0027 ±0.0027 



(28) 



IV. CONSEQUENCES FOR THE ISOVECTOR AND ISOSCALAR SPECTRAL FUNCTIONS 



If we drop terms second order in isospin breaking, then the p, u> and 
isoscalar vector current spectral functions are easily seen to be 



resonance contributions to the isovector and 
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[p 3 V)] y 

[p 8 V)]v 



F (3) 
P 



S P (s) 
''S u (a) 



p(8) 



(29) 



The results of the standard extractions, in contrast, are obtained by replacing Fp ' with F® M , F^ ' with y/iF^ M and 

F^ with \f3F^ M . The corrections to be applied to the standard contributions in order produce the true resonance 
contributions to the isovector and isoscalar spectral functions are then given by the ratios 

2 

= 0.979 ±0.0086 




= 1.189 ±0.065 



= 1.0054 ± 0.0054 . 



(30) 



From the above results we see that the standard procedure leads to an overestimate of the vector spectral function 
by 2.1 ± 0.9%. This is still noticeably smaller than the ~ 5% errors on the e + e~ — > hadrons cross-sections in the 
resonance region. As a result, it is not yet possible, when comparing to r data, to see the effect of the p 38 contributions 
to the e + e~ — ► hadrons spectral functions extracted using the conventional analysis. (See Fig. 10c of Ref. [1] for a 
comparison of the spectral functions as extracted from r and e + e _ data. The above correction would lower the e + e~ 
points by <~ 25 nanobarns in the region of the p peak.) A correction of this size, however, would certainly become 
important if one wished to make tests of CVC at the 1% level. 

A much greater surprise is the size of the correction required in the case of the u contribution to the isoscalar 
spectral function. While a 19% isospin breaking correction might sound unnaturally large, the size of the correction 
is, in fact, completely natural, and easily understood. The main features of the result follow from considering only 
the leading order contributions as discussed in section II above. Let us, therefore, consider the approximation in 
which one considers only the leading (0{m q , q )) contributions to p-uj mixing, neglects the "direct" isospin violating 
contributions to the vector meson decay constants (which are also higher order), and works in the ideal mixing/ SU(3) p 
approximation in which the EM decay constant of the pure isospin 1 component of the p is 3 times that of the pure 
isospin component of the u. Writing 



p = pi + e wi . 



u>i -epi 



(31) 



where the subscript I denotes the isospin pure states, and e is 0(6m), the physical EM decay constants become 



F 



p 

EM 



eFi 



I) 



3 
3e) 



(32) 



The fractional correction in the u> case is thus expected to be ~ 9 times as big as that for the p case. That the actual 
corrections turn out to be exactly a factor of 9 different is a numerical accident, but the large relative size of the 
corrections is completely natural, and associated with the smallness of the EM uj coupling and the pattern of mixing 
in the vector meson sector. Note also the fact that the corrections are of opposite signs is exactly what one expects 
based on the general arguments above. 

It is not just the isovector spectral function, with its relation to CVC, for which the corrections obtained above 
are of phenomenological interest. The difference of the isovector and isoscalar spectral functions also enters a number 
of interesting sum rules, and these sum rules must, therefore, also be corrected for the effects just discussed. As an 
example we will consider the extraction of the sixth (chiral) order low energy constant (LEC), Q{p), (in the notation 
of Rcfs. [29,9]) from the inverse moment chiral sum rule [29] 



r-OO 

j ds 

J 4ml 



(p™-p**)(s) 



16(m 2 K 



ml) 



3F 2 



Qip 2 ) 



1 



in-, 



+ 



Ll(^) + ^o(M 2 ) 
2tt 2 F 2 



ml log 



m- 



- m K log 



K 



(33) 
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In Eq. (33), L r k are the scale-dependent renormalized fourth order LEC's of Gasser and Letuwyler [17], and p is the 
ChPT renormalization scale. The form of this equation relies on the two-loop expressions for the 33 and 88 correlators 
obtained in Ref. [11]. The difference of the 33 and 88 spectral functions also enters a method of determining the strange 
current quark mass originally suggested by Narison [30]. In this application, a weighted integral over p 33 (q 2 ) — p 88 (q 2 ) 
is performed, the weight function being that which enters inclusive r decays. The corrections in this case, and the 
resulting values of the strange quark mass, will be treated in a separate paper [31]. 

In the analysis of the sum rule Eq. (33) performed in Ref. [29], the p contribution to the spectral integral was 
obtained from r decay data, and hence does not require the correction discussed above. The lo and <p contributions, 
however, are determined from the experimental V — ► e + e~ partial widths, and hence contain contributions from p 38 
which must be removed. The uncorrected p, lo and <j> contributions to the spectral integral are [29] 0.0374, —0.0103 
and —0.0204, respectively. Implementing the corrections above, one finds that the sum of these three contributions 
is reduced from 0.0067 to 0.0046 ± 0.0007, a downward shift of 31%. Including the estimates for the 4ir and KKtttt 
contributions as evaluated in Ref. [29], we find that the central value for Q(m 2 ) is shifted from 3.7 x 10~ 5 to 2.2 x 10~ 5 , 
a change of 41%. As noted above, because of the cancellations inherent in forming p 33 {q 2 ) — p 88 (q 2 ) (the combination 
vanishes in the SU(3)p limit), the effect of the isospin breaking corrections is large. A similar effect is found in the 
case of the strange quark mass analysis. 



V. SUMMARY 



We have shown that contributions to e + e~ — > hadrons involving an intermediate state p or intermediate state lo 
or <p contain contributions from the isospin violating 38 vector spectral function which arc not negligible, and must 
be removed if one wishes to extract the isovector 33 and isoscalar 88 spectral functions from e+e~ — > hadrons data. 
Using the results of a QCD sum rule analysis of the 38 correlator, we have been able to estimate the isospin violating 
vector meson decay constants required to make these subtractions. We find that the isovector spectral function is 
~ 2% smaller than what one would obtain by assuming it was identical to the full experimental p contribution, and 
that the lo contribution to the isoscalar spectral function is <~ 19% larger than what one would obtain from experiment 
without making this correction. We have also explained why it is unavoidable that (1) the isoscalar correction will 
be much larger than the isovector correction (by roughly an order of magnitude), and (2) the sign of the p and lo 
corrections in the isovector and isoscalar cases, respectively, will be opposite. A consequence of the second point is 
that all observables related to weighted integrals over the difference of the 33 and 88 spectral functions will receive 
large isospin breaking corrections, dominated by those which need to be made to correctly obtain the lo contribution 
to the 88 isoscalar term. 

Finally, we note that it might be possible to reduce the errors on the extractions of the isospin breaking decay 
constants by updating the sum rule analysis of the 38 correlator using recent improved values for the input parameters, 
and evaluating higher order a s corrections to the Wilson coefficients appearing in the D=2,4 terms of the OPE. This 
will be the subject of future investigations. 
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